Mapping spin-polarised transitions with atomic resolution 
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The coupling of Angstrom-sized electron probes with spin polarised electronic transitions shows 
that the inelastically scattered probe electron is in a mixed state containing electron vortices with 
non-zero orbital angular momentum. These electrons create an asymmetric intensity distribution 
in energy filtered diffraction patterns, giving access to maps of the magnetic moments with atomic 
resolution. A feasibility experiment shows evidence of the predicted effect. Potential applications 
are column-by-column maps of magnetic ordering, and the creation of Angstrom-sized free electrons 
with orbital angular momentum by inelastic scattering in a thin ferromagnetic foil. 
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I. INTRODUCTION 

With the availablity of electron vortices of sub-nm scale 
in the transmission electron microscope (TEM) [1-4] and 
its theoretical description [5Hl] many potential applica- 
tions come within reach, ranging from the transfer of 
angular momentum to nanoparticles, over utilisation of 
the intrinsic magnetic moment of vortex electrons to the 
probing of chirality[8j. Indeed, chiral electronic transi- 
tions were the first application of electron vortices in en- 
ergy loss magnetic chiral dichroism (EMCD) [Ij. 

The discovery of EMCD [9j was an unexpected alterna- 
tive to XMCD (X-ray magnetic circular dichroism) with 
the convenient side effect that additional information on 
the investigated material can be obtained simultaneously 
via standard analytical techniques [Tm [TT]. The spatial 
resolution of this electron microscopic technique is now 
in the nanometre range [12j. A modification of the tech- 
nique has been shown to be site selective [13j. The ex- 
cellent spatial resolution and the site selectivity are im- 
portant for the study of novel materials such as Heusler 
alloys [14] |T5], nanoparticles or interfaces [16]. Recent 
advances in electron microscopy have led to the imaging 
of condensed matter with subatomic resolution [T7HT9J. 
On this basis it has been speculated that the mapping 
of spin polarised electronic transitions - and thus the 
mapping of spin and orbital polarisation - on the atomic 
scale could be feasible in a TEM, applying the EMCD 
technique [20J. An incident plane wave affecting a spin po- 
larised L23 transition would break the mirror symmetry 
of non-magnetic transitions in the scattered wave. This 
symmetry breaking could be analysed either with asym- 
metric objective apertures or with a cylinder lens [2H [22 J. 
However, for technical reasons both approaches are un- 
realistic in the TEM. 

On second thoughts it becomes evident that the in- 
elastic interaction of an incident electron with a spin po- 



larised electronic transition creates a scattered electron 
with topological charge. One can thus apply the theory of 
vortex electrons [6j to the outgoing wave field. In a sense, 
this is a bottom-up application of the original idea of us- 
ing incident vortex electrons for EMCD. The reason that 
this works is the generalized reciprocity theorem [23j that 
confirms the equivalence of the incident and the outgoing 
(also called reciprocal [24J) electron for inelastic scatter- 
ing. This observation raises two questions: How can an 
EMCD signal be detected with Angstrom-sized Scanning 
TEM (STEM) probes; and can one produce electron vor- 
tices without holographic masks? 

Here, we present a theoretical and numerical analy- 
sis of the coupling between an Angstrom-sized STEM 
probe and an atom-sized vortex field via a chiral elec- 
tronic transition. It is shown that the inelastically scat- 
tered probe is in a mixed state containing electron vor- 
tices with non-zero angular momentum. These electrons 
create an asymmetric intensity distribution in energy fil- 
tered diffraction patterns, giving access to maps of the 
magnetic moments on an atomic column-by-column ba- 
sis. A feasibility experiment shows evidence of the pre- 
dicted effect. Finally potential applications are discussed: 
maps of magnetic ordering with atomic resolution; and 
the creation of free electrons with orbital angular mo- 
mentum and a diameter of about 0.1 nm by inelastic 
scattering in a thin ferromagnetic foil. 



II. THEORY 

We focus on the model of a thin (ideally one atom) 
layer of Fe. In this case the dynamical equation for the 
propagation of the probe's density matrix is considerably 
simplified. We give here only the basic equation for the 
propagator and refer the reader to the relevant litera- 
ture [20l[22l|25H2Il. 
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The inelastic intensity at energy loss E in the exit plane 
of the specimen is the diagonal term p(r, r) of the density 
matrix: 

p^(r,rO = J y'G,_,(r,x)GS_,(r^xOT^(x,x^z,zO 

where z' are variables along the optic axis and r, x 
are in planes perpendicular to the optic axis. G is the 
elastic propagator of electrons in the crystal. is the 
wave function of the incident electron at depth z and 
(lE = ko — ki is the minimum wave number transfer in 
the inelastic interaction, equal to the difference between 
wave numbers of the outgoing and the incident electron. 

The inelastic scattering kernel T£;(x, x.\ z, z^) is charac- 
teristic for the electronic transitions creating an energy 
loss E. For thin specimens and atomic columns with- 
out defects the z integration can be performed in closed 
form [28j . For spin polarised dipole transitions the atomic 
scattering kernel reads [22j 

tb(x,x)= yi e Cj'y^'- (2) 

/i=-l s=t,i 



with 

V±i(x) = 



27r Jo 



dq (3a) 



with the matrix element of the spherical Bessel function 

{jim = {iMQm 

between initial and final target states. = q'^-\-Qe with 
the characteristic wave number transfer = koE/ {2Eo). 
The coefficients (7^^ are weighting factors for spin-orbit 
coupling [S^^, "29"^ 30j, and n^'^ is the spin polarisation of 
the final state. The spin-orbit coupling of the initial state 
renders the coefficients C dependent on the total mag- 
netic quantum number j = l-\-s. For the L23 edges to be 
considered j = 1/2 or j = 3/2. The essential property 



that we will focus on is described by Eqj3a| It represents 
the outgoing inelastically scattered wave as an electron 
vortex with topology m = ±1 in the form of a Hankel 
transform that is easily computable from atomic wave 
functions. The azimuthal phase factor e^'^^ shows that 
the outgoing probe electron has orbital angular momen- 
tum. It should be mentioned that the probe beam is not 
spin polarised. It has acquired orbital angular momen- 
tum by spin-orbit coupling of the target electrons. 

The propagation of focussed probes through a thin 
specimen has regained interest in the context of real space 
STEM [31j. Even for elastic scattering the problem of 
propagating a focussed probe - as we shall adopt in the 



following - through a thin specimen to the detector poses 
considerable numerical problems. The inelastic interac- 
tion that can take place throughout the specimen adds 
another complexity. Therefore Eqjl] cannot be solved 
without approximations, at least with present numerical 
capacity. We restrict the discussion to a model system 
that allows to analyse the salient features of the inelas- 
tic coupling process with an accuracy comparable to the 
available experimental data. As such we choose a line 
of equally spaced atoms with given spin polarisation; we 
shall calculate the contributions from each transition in 
dipole approximation, discuss the signal from a single 
atom and finally build a line profile of the energy filtered 
signal from the array of atoms. For this model the elastic 
propagators G in Eqjl] collapse into delta functions and 
we have for the diagonal element (the density) 



PB(r,r)=T(r,r)(/.S(r)(/.o(r). 



(4) 



Eq. [2] allows us to disentangle the problem: Each transi- 
tion channel /i G [—1,1] can be factorised into a product 
of wave functions, and we get for the outgoing intensity 
in channel /i 



with 



and 



Pj ^{r,r) = Q^(^^(r)(^*(r) 



(5) 
(6) 

(7) 



The intensity in the diffraction plane is calculated from 
the 2D-Fourier transform of Eq. |5j 

Pjf,{ci,(i)=Cj^\FT,[^^{r)]\^ (8) 

as the trace over the 3 transition channels /i: 

/>j,£;(q,q) =Tr^[p^{ci,ci)] = ^/)^(q,q). (9) 

The trace operator shows formally that the scattered 
electron is in a mixed state. With the coefficients C for 
L23 transitions we can compute the outgoing intensity 
and the corresponding diffraction patterns. Without loss 
of generality we assume complete spin polarisation for 
the final target states (as is justified for the L edges of 
the 3d ferromagnets that we will use as a demonstration 
example). 

For complete spin polarisation the values for Cj^ are 
given m Table |l][52j. The image intensity Eqjl] will then 
contain different contributions from the scattering chan- 
nels for spin up and spin down polarisations. This differ- 
ence is the basis of EMCD. 

From here on we focus on the L3 edge and omit the 
index j for easier readability. Scrutinising Eq|5]for an in- 



cident plane wave, one notes that |(/?^| 



1^- 



P because 
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-1 1 


-1 1 


t 


0.056 0.111 0.167 
0.167 0.111 0.056 


0.278 0.222 0.167 
0.167 0.222 0.278 


unpolarised 


0.111 0.111 0.111 


0.222 0.222 0.222 



TABLE I. Prefactors C]'^ . The first two rows are the weigth- 
ing factors for the transitions when the final states are com- 
pletely (up or down) spin polarised. The third row gives the 
weighting factors for unpolarised final states. (All per elec- 
tron) . 



the phase factors of both the incident plane wave and the 
kernel - e^*" - cancel in the intensity. That means that 
there will be no difference in the image predicted by Eq|5] 
for spin up or for spin down polarisations. The same is 
true for the intensity in the diffraction plane Eq|9] Es- 
sentially, it is not possible to see spin polarisation from 
single atoms in the TEM without further action. (In 
the standard EMCD geometry, one uses the interference 
terms caused by Bragg scattering of the outgoing atomic 
vortices on the lattice, which is different for spin up and 
spin down polarisations.) 

The situation changes when using a STEM probe in- 
stead of a plane incident wave. We analyse the situ- 
ation qualitatively before considering numerical simula- 
tions. For the explanation of the effect we assume a nar- 
row focussed probe given by the Airy function A with 
a diameter much smaller than the distance R from the 
atom, essentially so small that that the amplitude of the 
kernel is almost constant within the probe. (In the nu- 
merical simulation this condition is relieved.) Then we 
can approximate the outgoing wave in channel ji as 

^^(r) = tA^(r - RM(r) « ^(r)|^^(i?)|e^''"W. (10) 

Note that the outgoing wave is more extended than the 
Airy disk because of the long range Coulomb coupling 
force; in Fig. [l] it is drawn not to scale (even larger for 
better visibility). The azimuth angle within the outgo- 
ing disk is a = a(R) + l-iy/R in the coordinate system 
shown in Fig {T] and we see that the phase of the outgoing 
wave depends on the position of the STEM probe and 
changes sign when going from /i = 1 to -1. Via the shift 
theorem the diffraction pattern will be proportional to 
the Fourier transform of the Airy function shifted in qy 
direction by — 1/i?, 0, 1/i? for the 3 transition channels, 
that are shifted disk functions II (q) 

= ^jMlV^M(^)l'n(q + /iq^/i^). (ii) 

Having established an observable that has the signature 
of a particular transition channel the spin polarisation of 
a single atomic column can be determined. 



III. NUMERICAL SIMULATIONS 

Simulations were performed for an incident probe of 
100 keV. First we construct the outgoing signal for each of 
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FIG. 1. Schematic of the principle of probe- vortex coupling, 
a) Top view of a narrow incident Airy disk (gray), focussed at 
a distance R in scan direction (x) from the atom (red). The 
scattering kernel Te is symbolised as a diffuse cloud with a 
color coded phase (hereafter called rainbow wheel), increasig 
from blue to red in clockwise rotation for the ii — 1 channel. 
According to EqjlO] the outgoing wave (gray circle) has ac- 
quired a phase ramp da = liy/R in the chiral trasition. b) 
The phase ramp translates into a shift of the diffraction disk 
by 11 /R in direction qy. c) the qx extension is squeezed into 
one pixel on the detector in the (q, E) geometry. Note that 
the energy loss axis must be perpendicular to qy , which is pro- 
portional to the scattering angle, qy — koO. The (0, E) map 
is a true image taken on magnetite ranging from the oxygen 
K edge at - 540 eV to the Fe L23 edge at - 710 eV. 



the 3 transition channels, p^(r, r), shown in Figj2j These 
are L3 energy filtered images of a STEM probe of 0.1 nm 
diameter scanning across a single atom. From top to bot- 
tom are the transition channels ji = —1, 0, +1. From left 
to right the distance R to the atom is -2,-1,0,1,2 at. units. 
Note that the STEM probe is always in the center of the 
images. The side structures at distances of = ±2 at.u. 
are signals from the second maximum of the Airy disk 
that coincides here with the atom centre. Brightness 
codes intensity of the image, color codes for the phase 
of the wave function. (Color online: rainbow colors from 
— TT to tt). When the Airy disk sits on the atom the outgo- 
ing beam is a true atomic vortex with topological charge 
/i G [— 1,0, 1]. At larger distance phase ramps in the 
Airy disks develop, visible as continuous color variations. 
They change sign when crossing the atom centre and are 
opposite for /i = ±1. Each square has a side length of 5 
at. units (0.26 nm). 

These phase ramps are responsible for the correspond- 
ing shifts of the diffraction patterns shown in FigjSj The 
shift of the patterns in vertical direction (qy) is opposite 
for /i = ±1 and depends on the position of the probe. 
The /i = channel does not show any shift because 
it lacks a phase ramp in the image. The patterns are 
smeared by convolution of the incident probe disk with 




FIG. 2. Fe L3 energy filtered real space exit wave functions 
Lp^ — Eq. [7| — of a 100 keV STEM probe of 0.1 nm diam- 
eter scanning across a single atom for the 3 dipole allowed 
transition channels. From top to bottom: transition channels 
jjL = —1, 0, +1. From left to right: distance R to the atom -2,- 
1,0,1,2 at. units. Brightness codes for intensity of the image, 
color codes for the phase of the wave function (Color online: 
rainbow wheel as given in Fig.[T]). When the Airy disk sits on 
the atom (center column) the outgoing beam is a real vortex 
with /i G [—1,0,1]. At a distance phase ramps in the Airy 
disks develop, changing sign with R and with /i. Each square 
has a side length of 5 at. units (^ 0.26 nm). 
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FIG. 3. Fe L3 energy filtered diffraction patterns p(q, q)/^ 
— Eq. [8] — of a 100 keV STEM probe of 0.1 nm diameter 
scanning across a single atom, corresponding to Fig [2] From 
top to bottom: transition channels /x = — 1,0,+1. From left 
to right: distance R to the atom -2,-1,0,1,2 at. units. Note 
the shift of the patterns for /x = ±1 in vertical direction qy as 
predicted in Figjl}), depending on the position of the probe, 
and the inversion of shifts with change of sign. The fi — 
channel does not show any shift because there is no phase 
ramp in the image. Convergence angle 18 mrad. The images 
have a side length of ±50 mrad. 



the inelastic scattering kernel that has an extension of 
SO - AE/{2Eo) ^ 3.5 mrad. 

Since the /i = ±1 channels contribute differently for 
spin up/spin down polarisation via the coefficients Cj^, 
Eqj6j the diffraction patterns will be different for these 
two cases. (Note that this is not the case in the fil- 
tered real space image: there, the difference is only in 
the phase, not in the intensity distribution). 

Monitoring energy filtered diffraction patterns of a 
scanned probe means measuring a multidimensional data 
cube because background subtraction and multiple scat- 
tering deconvolution of energy loss spectra require a 
range of losses. Such data has two dimensions (qx,qy) 
in the reciprocal space, one in the energy-loss (eV), and 
one (x) or two (x,y) in the real space, depending on the 
scanning pattern. This creates huge data files - a scan 
over 1 elementary cell in magnetite with 0.02 nm step 
width with 256^ pixels in the diffraction pattern would 
give ~ 60 Mb per energy channel, approaching ~ 10 GB 
for a whole spectrum - and is impractical. One can how- 
ever exploit a remarkable feature apparent in the simu- 
lations: The diffraction patterns show only asymmetry 
with respect to the coordinate g^, which is the Fourier 
transformed variable of y. An experimental setup could 
then discard the Qx variable by integration without in- 
formation loss. This is exactly what is realised in the 
(g, E) geometry. There, the Qx axis is squeezed onto one 
pixel of the specrometre by compressive lenses whereas 
the Qy axis is retained, being projected on the detector 
perpendicular to the energy loss axis. This is sketched in 
Figjlp which is a (qy^E) data set for a fixed positon R 



. Selecting the L3 white line from the whole (qy^E^R) 
data, a reduced subset with axes (i?, qy) is obtained. The 
scan direction must be perpendicular to y in the coordi- 
nate system of Fig|l] Contrary to XMCD or EMCD, the 
new technique operates on a single white line only (the 
stronger I/3 line here). This is important because often 
the L2 edge is too faint to obtain sensible results. (It 
should be noted that for the separation of spin and or- 
bital moments both L2 and L3 edges are needed [32 , 33]. 

These subsets can be constructed from the previous re- 
sults, simply by integrating the intensity of Figj3]over the 
"'squeezed'" variable qx and summing over the 3 transi- 
tion channels in Eqj9j The result is shown in Figj4] for 
spin up and spin down configurations. The asymmetry of 
the intensity distribution with respect to the atom posi- 
tion at = is indicative for spin polarised transitions. 

To enhance the asymmetry Figj4] suggests to avoid 
the central part around qy = which adds only spin- 
insensitive intensity thus increasing the noise level. In- 
tegration of the density matrix EqjSjover the scattering 
angle in the top and bottom parts yields two scans over 
a single atom 

/pqy2 
dqx / p(q,q)c^<72/ 

/r-qyi 
dqx / p{q,q)dqy (12) 
-qy2 

where qy = ko6. This is shown in Figji]: for 6yi = 
20 mrad, 0y2 = 50 mrad. The difference in position be- 
tween the maxima in the two scans is ^ 0.06 nm, indi- 
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FIG. 4. Fe L3 energy filtered {R,qy) signal distribution of 
a STEM probe of 0.1 nm diameter scanning across a single 
atom, corresponing to Fig|2] a: spin up polarisation of the 
atom; b: spin down polarisation. The top/bottom asymmetry 
is well visible and can be used to determine the spin polarisa- 
tion, c: Full line: Intensity integrated over the upper detector 
half (20 to 50 mrad) hne scan of a STEM probe of 0.1 nm di- 
ameter across a single atom, corresponding to a). Dashed 
line: same for the lower detector half (-50 to -20 mrad). 



eating that Cs corrected machines and extreme stability 
are needed to see the effect. Even so, the signal will 
be very faint, such that noise will tend to override the 
effect. Dynamical diffraction of the incident and the out- 
going electron on the lattice and remaining aberrations 
of the probe forming lens and the spectrometer will also 
complicate the situation. Quantitative spin detection will 
therefore need elaborate calculations. 



IV. EXPERIMENTAL EVIDENCE 

We have performed a feasibility experiment at the 
Daresbury SuperSTEM facility using a NION Ultra- 
STEMIOO microscope [34J. The instrument is a dedi- 
cated, aberration corrected STEM operated at 100 kV 
with a cold-FEG emitter. Its 3rd generation C3/C5 aber- 
ration corrector allows a typical probe size of 0.1 nm at 
a beam current of 30 pA. EELS spectrum imaging at 
atomic resolution is done with a DigiScan2 scanning unit 
in combination with a Gatan Enfina spectrometer. The 




FIG. 5. HAADF image of the border of the magnetite 
nanoparticle in [111] zone axis orientation. The scan line is 
indicated as a box, the axis of energy dispersion (arrow) is al- 
most perfectly parallel to the scan direction, thus qy\\y. The 
bright dots are the A columns containing 3 Fe atoms per unit 
cell in this projection. Their projected distance is 0.33 nm. 



difficulty of the experiment stems from the extreme de- 
mand on stability, the rather long dwell times causing 
beam damage, remaining aberrations, and the low count 
rates. Spectra were collected during 0.2 s per position, 
corresponding to a dose of ~ 5. 10^ electrons focussed in 
the STEM spot. 

We investigated a platelet-like magnetite nanoparticle 
of ~ 15 nm diameter, FigjS] The thickness in this re- 
gion was between 5 and 10 nm, resembling as close as 
possible the single atomic row model. The STEM probe 
was scanned over a line of atoms marked in figjS] as a 
rectangle. Magnetite is an inverted cubic spinel XY2O4 
with Fe^^ ions at tetrahedral X sites. The octahedral Y 
sites are randomly occupied by Fe^+ and Fe^^ ions. X 
and Y sites are antiferromagnetically coupled. In (111) 
zone axis there are two types of Fe columns (hereafter 
called A and B). A contains two X sites and one Y site 
per elementary cell, B contains only one Y site. The B 
columns are not visible in the HAADF image FigjS] be- 
cause of lower scattering strength and dynamical diffrac- 
tion, confirmed by multislice simulations. However, very 
faint side maxima from the B columns can be seen in 
Fig. [6^. Shown there are the top and bottom profiles 
of the scan integrated from 20 to 50 mrad after stan- 
dard background subtraction, drift correction and usual 
removal of the continuum signal beneath the white line. 
The post-edge continuum, often used in standard EMCD 
for normalising, did not show periodic variations. There- 
fore it was not necessary to correct for, the more so as 
this would have induced additional noise. A shift of the 
traces with respect to the atomic positions can already 
be guessed although the noise is almost overriding the 
signal. The number of electrons collected in the I/3 edge 
was ~ 250 in each half detector, causing a theoretical 
shot noise level of 3a = 47. Pre-edge extrapolation, in- 
stability and other error sources add noise such that the 
3<j relative error amounts to ~ 20 — 30%. 

In order to demonstrate the predicted effect qualita- 
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FIG. 6. a: Scan of the Fe L3 white hne signal after standard 
background subtraction, drift correction and removal of the 
continuum. The filled curve shows the scan using the upper 
detector half (20 to 50 mrad), the empty curve is a scan using 
the lower detector half (-50 to -20 mrad). The HAADF signal 
is superposed (dashed) with maxima indicating the positions 
of the A columns. The B columns are visible as faint sub- 
sidiary peaks, b: Same scans after Fourier filtering. Atom 
positions (MAADF maxima) are marked with vertical lines. 

tively a Fourier analysis was perforraed on the two scans, 
retaining only coefficients up to lattice periodicity. The 
result is shown in Fig.|6]3. As predicted in Figj4]the max- 
ima are shifted to both sides of the atom centres that are 
marked by vertical lines. The average distance between 
the top and bottom scan maxima is 0.076 nm. One ob- 
serves a rather large shift at the 0.33 nm position that 
is caused by some irregularity (probably sudden drift). 
Excluding this value, the average difference between the 
respective maxima in the two scans is 0.65 nm, in good 
agreement with the simulational result of 0.6 nm seen in 
Figj4]3. Despite the crude approximations and the simple 
model the agreement is surprisingly good. More accurate 
models can be devised, but this is beyond the scope of 
the present paper as noise, drift and beam damage pose 
narrow limits on the interpretation of the data. More 
elaborate experiments with ultrathin magnetic specimens 
must be performed in order to confirm the present find- 
ings. 



V. POTENTIAL APPLICATIONS 

A. Detection of magnetic order 

As described above and shown in Fig|4j the 
top/bottom asymmetry in the (R^Qy) data set is caused 
by the spin polarisation on the atomic site. That pro- 
vides a method for spin mapping. A straightforward way 
to do so is to take the first derivative of the difference 
(top-bottom) line scans (p^ — p'_). Since the slope of 
the difference signal is strongest at the atomic sites, this 
scan gives directly the sign and position of the magnetic 
moment. Here we assume a hypothetical system with the 
same lattice constant as magnetite, also in (111) zone axis 
orientation, to make connection to real systems. FigjT] 
shows the derivative {pj^ — p-)' {R) along a line scan as 
in FigjSj for different assumptions of the magnetic or- 
dering (ferromagnetic, antiferromagnetic, and ferrimag- 
netic). To take account of the dechanneling and defo- 
cussing throughout the specimen the data was convolved 
with a 0.1 nm broadening function. The last panel re- 
sembles magnetite, in fact. In (111) zone axis projection 
the A columns contain 2 atoms spin up and on atom 
spin down, and the B columns contain one atom spin 
down, oriented along the z axis in the magnetic field of 
the objective lens. On the other hand, the channeling 
is stronger on the deeper potential (A columns contain- 
ing 3 atoms as compared to 1 on B) so the beam will 
see more from the A columns than from the B columns. 
This will weigh the A columns stronger than the B ones. 
The exact weighting factor is impossible to obtain with- 
out solving the dynamical equation for the propagator 
in the lattice, but a weighting coefficient between -50% 
and -25% for the B columns is reasonable. Here, -30% 
were assumed for the simulation. Theoretically, it should 
be possible to perform such scans not only along lines 
but over areas. This should give atom-resolved maps of 
the element specific magnetic moments and the magnetic 
ordering. But before establishing such an analytical tech- 
nique the problems related to noise, stability, aberrations 
and dynamical diffraction must be solved. 

B. Free electrons with angular momentum 

If an Angstrom sized spot is focussed exactly at the 
centre of an atom in a ferromagnet the scattered electron 
has acquired orbital momentum. It is important to note 
that it is in a mixed state with contributions from the 
three transition channels p G [—1,1], each creating a pure 
vortex state with topological charge p. The expectation 
value of the angular momentum can be calculated from 
the coefficients Cj^ given in Eq. |6] 

(^.) = ^^^^^f^ G[-0.167n,0.167n] 

depending on the spin polarisation of the atom. This 
is a unique and simple method to create free electrons 
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FIG. 7. Hypothetical line scans showing the first derivative 
d{p-^{R) — p-{R)/dr over an atomic row as in Fig. [s] for a: 
ferromagnetic, b: antiferromagnetic and c: ferrimagnetic or- 
dering, c) was calculated with a relative strength of -30 % for 
the spin down moments, thus taking into account the strong 
channeling at the A columns. 




FIG. 8. The phase and amplitude of the incident Airy disk 
of 0.1 nm diameter (left) and for the scattered wave from the 
fi — 1 transition. The atom is displaced from the beam center 
to the left by 0, 0.01, and 0.05 nm. The images have a side 
length of 0.2 nm. 




FIG. 9. Intensity of the incident Airy disk (left) and Tr[p] for 
the scattered wave (mixed state from 3 transition channels) 
for atom displacements as in Fig. [8] 



beam for the /i = 1 transition when the atom is in the 
center, 0.01 nm, and 0.05 nm sideways. The vortex struc- 
ture disappears rapidly, also visible in the intensity dis- 
tribution Fig. [9] Here, the total electron density Tr^[p] is 
shown. The central dip, characteristic for the topological 
charge, remains at a deviation of 0.01 nm but has disap- 
peared for 0.05 nm. Fig.[T0^ compares the incident Airy 
disk with the radial profile of the outgoing mixed state. 
The scattered state is even narrower than the Airy disk. 
The high sensitivity of the outgoing vortex state to the 
probe position could be used for a more direct method of 
spin mapping with sub-atomic resolution than described 
above, e.g. with a vortex filter such as a holographic 
mask. 

We simulated also the scattered state when the inci- 
dent beam is broader (0.5 nm diameter). The result is 
shown in Fig, [ToJd. The outgoing beam is broader than 
in case a) but still much narrower than the incident Airy 
disk, and is again a superposition of electron vortices. 
The not so surprising consequence is that a thin ferro- 
magnetic foil in an electron beam creates free electrons 
carrying angular momentum after energy filtering. This 
shows that vortices always have been there in EELS ex- 
periments on magnetic materials. 



VI. CONCLUSION 



with orbital momentum, although the efficiency and the 
available momentum is probably too low to manipulate 
nanoparticles or even atoms by the torque. Another 
problem is the precise positioning of the probe on the 
atom. Fig. |8] shows the phase and amplitude of the inci- 
dent beam of one Angstrom diameter and the scattered 



The coupling of an Angstrom-sized electron probe to 
a spin polarised transition creates a mixed state that 
contains electron vortices with non-zero orbital momen- 
tum. These states break the symmetry of the scattering 
distribution in the far field in a way characteristic for 
the chirality of the transition, a fact that can be used 
for the imaging of electron spins in real space with sub- 
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FIG. 10. a) Radial profiles of the incident one- Angstrom 
probe (dashed) and the outgoing electron density of the mixed 
state when the atom is exactly centred on the beam, b) The 
same for an incident probe of 5 Angstrom diameter. 



Angstrom resolution. A tentative experiment on a mag- 
netite nanoparticle shows the expected asymmetry. 



Apart from probing the local magnetic ordering, im- 
portant for a number of technologically promising mate- 
rials such as Heusler alloys, the proposed method bears 
promise for the mapping of spin polarisations of single 
atomic columns, be that in the vicinity of interfaces, mag- 
netically dead layers, or magnetic core-shell structures. 



The creation of free electrons carrying angular momen- 
tum is theoretically feasible via spin polarised electronic 
transitions. This works even for relatively broad incident 
beams passing a thin ferromagnetic foil. 
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